Astronomy & Astrophysics manuscript no. 
(will be inserted by hand later) 



Dynamical modeling of large scale asymmetries in the (3 Pictoris 
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Abstract. We report a new and complete model of the /3 Pictoris disk, which succeeds in accounting for both the 
surface brightness distribution, warp characteristics, the outer "butterfly" asymmetry as observed by HST/STIS 
in scattered light, as well as the infrared emission. Our model includes the presence of a disk of planetesimals 
extending out to 120-150 AU, p erturbed gravitationally by a giant planet on an inclined orbit, following the ap- 
proach of Mouillet et al. (1997b). At any time, the planetesimal disk is assumed to be the source of a distribution 



of grains produced through collisional evolution, with the same initial orbital parameter distribution. The steady 
state spatial grain distribution is found incorporating the effects of radiation pressure which can cause the distri- 
bution of the smallest particles to become very distended. With realistic assumptions about the grains' chemical 
properties, the modeling confirms the previously evident need for an additional population of hot grains close to 
the star, to account for the 12 fim fluxes at short distances from the star. It also indicates that this population 
cannot explain the outer 12 /im flux distribution when the effects of gravity and radiation pressure determine the 
distribution. Very small grains, produced by collisions among aggregates, are tentatively proposed to account for 
this 12 /xm outer emission. 

Key words. Stars: circumstellar matter - Stars: /3 Pictoris 



1. Introduction 

The (3 Pictoris gaseous and dusty disk has been exten- 
sively studied for 15 years. The dust disk, seen nearly 
edge-on, extends to at least a distance of 100 AU from the 



are best fitted by the emission of small (submicronic) crys- 
talline silicates (Knacke ct al. 1993; Li & Greenberg 1998). 
It was suggested that these grains could result from comet 
evaporation, a scenario that had already been proposed to 



central] star, with a sharp decrease in the surface brigliL- 



explain the very peculiar spectroscopic activity (Lagrange 



ness distribution of scattered tight beyond about 120 AU. 



The dynamics of the small grains producing the scattered 

light is determined by radiation pressure, collisions with 
other grains and/or evaporation and to a lesser extent, 
Poynting-Robertson drag. The images at 12 /im show that 
the maximum of the dust surface density distribution is 
located between 80 AU and 100 AU flPantin et al. 1997|, 



et al. 1987; Ferlet et al. 1987). Furthermore, Lecavelier de; 
Etangs et al. (1996| ) proposed that the appearance of the 



entire disk, as seen in scattered light, could be explained 
by grains released from evaporating comets at distances 
of 15-30 AU and subsequently pushed out by radiation 
pressure. 



assuming th e Hipparcos distance of 19.28 pc from |Crifo| 1.1. Asymmetries in the dust disk 



et al. (1997) for (3 Pictoris). The total dust mass, as mea- 
sured from sub- millimeter data, ranges between a few and 
a few tens of lunar masses (Zuckcrman and Bccklin 1993; 

20001) ~ 



Holland et al. 1998; Dent et al 



The 10 /im silicate spectrum flKnacke ct al. 1993| ) ap- 
pears to mimic that of Halley's comet but not that of 



comet Hale-Bopp (Lagage et al. 199E). The 10 /im features 
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The Northeast and Southwest extensions of the dust disk 
have been found to be asymmetric in scattered light as 
well as in thermal emission. At visual wavelength, Kalag 
|fc Jewitt (1995 ) noticed five radial or vertical asymme- 
tries beyond 150 AU. In the Northeast extension for ex- 
ample, the vertical brightness distribution in scattered 
light is asymmetric with respect to the disk mid-plane 
out to hundreds of AU and this situation is reversed in 



the Southwest extension ("butterfly" asymmetry). Kalas 
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et al. (2000) and Larwood et al. (2001) recently proposed 
that this asymmetry is produced by a close encounter with 
a MOV star in the last 10 5 years. 

In the inner part of the disk (at about 60-70 AU) the 
12 /im images reveal a factor of 3 side to side brightness 
asymmetry (Pantin et al. 1997) which is much smaller 



at optical and near- infrared wavelengths (Mouillet et al. 
1997aQat most a factor of 1.5). Orbiting evaporating bod- 
ies on eccentric orbits with the same longitude of perias- 



tron can produce such axial asymmetries (Lecavelier des 



Etangs et al. 1996; Lecavelier Des Etangs 1998). 

Adapti ve optics ([Mouillet et al. 1997b|) and HST ob- 
servations (Burrows et al. 1995; Heap et al. 2000) revealed 
the inner warping of its dust disk at about 70 AU and in 



the same direction as the "butterfly" asymmetry. Mouillet 



et al. (1997b) modeled this warp as being due to grains 
which followed the distribution of a parent body plan- 
etesimal disk which was perturbed by the gravitational 
influence of a Jupiter-like planet on an orbit inclined at 3 
degrees to the disk mid-plane. 

1.2. Propagation of the warp 

The important effect of perturbation by the giant planet 
is to cause the parent body disk to precess differentially. 
The more rapid precession in the inner regions causes 
the parent body disk to become coplanar with the plane- 
tary orbit. At large distances the inclination of the parent 
body disk retains its initial value. The warp is located 
in the transition region where the local orbital precession 
timescale is comparable to the age of the system. The lo- 
cation of the warp accordingly propagates outwards with 
time. 

Mouillet et al. (1997b) found an approximate relation- 
ship between the position of the warp, the mass ratio of 
planet to star, its orbital (assumed circular) radius and 
the age of the system in the form: 



t 



0.2 



Here M», M, R w , D, t and t un n are the mass of the star, 
the mass of the planet, the radius of the warp in the parent 
body disk, the planetary orbital radius, the age of the sys- 
tem and the time unit \J (10AU) 3 / (GM*) ~ 5.2 y respec- 
tively. For example, to obtain R w = 70AU, M /M* = 10~ 3 
and D = 10AU requires the age of the system to be 
~ 2.1 x 10 7 y. Such orbital parameters for the proposed 
perturbing planet allow the spectroscopic activi ty to be 
successfully reproduced ( Bcust fc Morbidelli 2000 , and ref- 
erences therein). 

Note that when, as assumed here, the warp location is 
significantly larger than the orbital radius, that location 
depends only on the product of the age of the system, 
the square of the orbital radius and the planetary mass 
ratio. We comment that, because it does not incorporate 
a non axisymmctric perturbation, the model proposed by 
Lecavelier des Etangs et al. (1996| ) cannot account for the 



inner warp, but the influence of radiation pressure, clearly 
revealed by their study, is certainly of great importance. 

1.3. Aims of the present modeling 

Complete modeling of the disk has not been performed so 
far. In this paper we formulate a model of the disk which 
reproduces both the scattered light images, in particular 
the warped disk and the "butterfly" asymmetry, and also 
the infrared (IR) images and fluxes. In Section ||, we out- 
line the basis of this model and Section || presents the 
scattered light images which can be compared with the 
recent HST/STIS data ( [Heap et al. 2000| ). In section f§ 
we synthesize the IR 12 fim image and compare with the 
observations. In section |5| we summarize and discuss the 
proposed model in a more general context. 



2. The basic model 

2.1. The parent body disk 

The short lifetime of the grains, in comparison to the age 
of the system, as a result of the action of collisions and ra- 
diation pressure leads to the hypothesis of the existence of 
a parent body disk consisting of planetesimals with sizes 
larger than 1 m which acts as a source for the grains as col- 
lision products (see for example the discussion in Mouillet 
|et al. 1997b , and references therein). 

As in Mouillet ct al. (1997b| ) we suppose that the phys- 
ical collision time in the parent body disk is sufficiently 
long compared to the age of the system that it can be ap- 
proximated as a collisionless disk lying between 1.5D and 
15D. The disk aspect ratio H/r is taken to be 0.1. Through 
our modeling procedure we are able to construct arbitrary 
surface density profiles. This is done by constructing par- 
ent body disks consisting of weighted contributions from 
100 elementary rings composed of 10 4 particles and with 
a ratio of inner and outer radii equal to 10 1 /" ~ 1.0235. 
This is possible because the particles are collisionless. 



2.2. Orbital parameters of the perturbing planet 

For working parameters we adopt D — 10AU, being also 
the model unit of distance, planet to star mass ratio 
M/M* = 10~ 3 , planetary orbital eccentricity of zero with 
an inclination of 3° with respect to the initial mid-plane 
of the parent body disk. However, as above, we note that 
the important parameter is MD 2 rather than M and D 



separately. Mouillet et al. (1997b) also show that the prop- 



agation of the warp is unaffected for eccentricities of the 
planetary orbit up to 0.5. however, large eccentricities may 
lead to side to side (radial) asymmetries of the parent 
body disc which in turn may lead to some asymmetries in 
the generation of small dust particles. The study of asym- 
metries produced in this way is beyond the scope of this 
paper. 
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2.3. Parent body disk precession 

The parent body disks are modeled here using 10 6 parti- 
cles. It is not possible to carry out orbital integrations over 
the system lifetime for such a number of particles. Instead, 
after realizing the parent body disk by allocating particles 
a random inclination and eccentricity consistent with the 
aspect ratio, the system was allowed to phase mix by inte- 
grating for 20 outer orbital periods. At this stage the or- 
bital plane of each particle was precessed at the local rate 
for the age of the system while each particle retained its 
orbital phase. Finally the system was phase mixed again. 
This procedure is equivalent to evolving the orbits using 
the time averaged potential of the planet and its validity 
has been checked by comparison with complete integra- 
tions over time spans where these were feasible. The local 
precession rate for near circular orbits was taken to be 



( pviouillet et al. 1997b| ) 
3GMD 2 



Here r is the orbital radius and = -JGM^/r 3 is the lo- 
cal Keplerian frequency. The age of the system adopted 
here, corresponding to the mass ratio M/M* = 10~ 3 , 
was 2.6 x 10 7 y co nsistent with the recent estimated age 
of /3Pictoris from Barrado y Navascues et al. (1999 ) and 
which yields R w ~ 75 AU. But note as above that the sig- 
nificant parameter is the product of the age and the mass 
ratio. 

The configuration of the warped parent body disk, ob- 
tained using the above procedures, with the surface den- 
sity we adopted is illustrated in Figure |[ 



due to radiation pressure to that due to the gravity of the 
central star, /3 pr = F ra d /F glav . The case pl = corre- 
sponds to a bound dust particle with the same motion as 
its PB. When 3 pr = 0.5, particles have zero energy and 
so are just unbound. We consider 11 values of /3 pr equally 
spaced between and including 0.45 and 0.001. For each 
value of /3 pr , a dust distribution is generated from the 
PB by phase mixing 10 6 particles for 20 orbital periods 
measured at the outer radius of the disk. The time-scales 
required for phase mixing are always short compared to 
the time-scale for propagating the warp. Also the optical 
depth of the disk is large enough that the timescale for 
destruction of the dust by collisions is shorter than any 
dra g timescale due to the Poynting-Robertson effect (see 
also Mouillct ct al. 1997b|) . The Poynting-Robertson effect 
then is neglected here. 



2.5. Modeling of the scattered light and thermal 
images 

2.5.1. Efficient modeling of the PB surface density 
distribution 

To construct the final image of the disk in scattered light 
or in thermal emission arising from any PB surface density 
distribution in practice, a radial weight function is applied 
to the contribution from each dust particle depending on 
the elementary ring between 1.5D and 15D from which 
it originated (the distance to this ring also corresponds 
to the periastron distance of the dust particle). Several 
S(r) distributions can therefore be investigated using this 
technique in a reasonable amount of time. 




Fig. 1. Vertical position of the mean PB disk mass as- 
suming the surface density distribution of Figure |. The 
vertical and radial scales are not the same in this Figure. 
The warped disk extends out to 140-150 AU with a maxi- 
mum close to 75 AU from the star and a vertical extension 
of ^4 AU to the initial mid-plane of the PB disk. 



2.5.2. Grain properties 

The image of the disk also depends on the physical prop- 
erties of the grains. Here we have adopted the comet dust 
model proposed by Li fc Greenberg (1998 , and references 
therein) to reproduce the (3 Pictoris spectral energy distri- 
bution (SED). In this model, the dust consists of porous 
aggregates made fro m a silicate core co ated with organic 
refractories (see also Pantin et al. 1997 ). A mor e detailed 
account of the grain model used here is given in Augereai; 
|et al. (1999 ). Initially, we fix the porosity of the grains 



to be 0.95 and assume that most are amorphous (Li & 
iGreenberg 19981 ; [Lagage et al. 1999[ ). 



2.4. Simulation of the dust disk 

We assume that each planetesimal (parent body or here- 
after PB) is the progenitor of small particles (produced 
for instance by collisions) that are significantly affected 
by the radiation pressure of the central star. We assume 
no initial relative velocity between the PB and the small- 
est particles produced. A consequence of this assumption 
and the near circular motion of the PB is that the perias- 
tron of a dust particle equals that of the PB. The effect 
of radiation pressure is quantified by the ratio of the force 



2.5.3. Presence of ice 

Water ice may be present on the surfaces of the ele- 
mentary submicronic particles that compose an aggregate 
( Greenberg and Hage 199C| , and references therein). The 
amount of ice is quantified by a volume percentage in ad- 
dition to that of vacuum (due to the porosity). A grain is 
assumed icy if it is produced further than the ice subli- 
mation distance. That distance depends on the grain size 
and amount of ice (typically between 20 AU for grains 
larger than ~ 100 /im and 100 AU for grains smaller than 
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~ 1 /mi). Due to the effect of radiation pressure, non-icy 
grains can be present at very large distances as well as icy 
grains, possibly mixing two dust chemical compositions at 
a given distance from the star. 



Because dust particles trace the distribution of the PB 
disk (Mouillct et al. 1997a), qualitatively, these remarks 
imply that most of the PB disk mass is contained within 
an annulus lying between ~50 AU and 120 AU. 



2.5.4. Grain sizes 



3.2. Surface brightness distribution 



For the range of /3 pr values we consider, there is a one-to- 
one relationship between the value of /3 pr and the grain 
s ize a which can b e expressed in the form a = Kj3~^ 
( Artymowicz 1988| , for instance), where K depends on the 
optical properties of the grain which arc computed using 
the Mie theory Q ( [Bohrcn fc Huffman 1983) ). borne exam- 
ples of K values for different amounts of ice and porosities 
can be derived from the values of a^j pr= o.5 given in Table 
|l|. A grain size distribution with dn(a)/da oc a~ 3 ' 5 is as- 
sumed to take into account the relative number n(a) of 
grains with radius greater than a produced in the disk by 
collision processes. For each /3 pr value (or grain size), the 
appropriate relative number density and thus the contri- 
bution to a scattered light or thermal emission image is 
computed. The 11 independent images then are correctly 
superimposed with appropriate weight factors so as to pro- 
duce a final image of the disk for any required PB surface 
density distribution or wavelength of interest. 



2.5.5. Orientation of the disk to the line of sight 

In order to complete the above procedure, the disk orienta- 
tion with respect to the observer has to be specified. This 
introduces two additional parameters: the line of nodes of 
the perturbing planet and the disk inclination with respect 
to a plane normal to the line of sight. 



3. Scattered light images 

3.1. PB surface density distribution 

The first stage in the modeling procedure is to find a sur- 
face density distribution for the PB disk that leads to 
consistent results with both scattered light observations 
and thermal emission data. The range of possible surface 
density distributions can be restricted using the following 
observational constraints: 



1/ Both the resolved 12/tm images (Pantin et al. 1997) 
and SED modeling (Li & Greenberg 1998) indicate that 
little dust lies inside ~50 AU. 

2/ Beyond about 120 AU, the radial surface brightness 
profile in scattered light revealed by STIS is roughly pro- 
portional to r~ 5 (Heap et al. 200C). For an edge-on disk 



like the (3 Pictoris one, this is in very good agreement with 
the profile expected from small dust particles produced by 
colliding large bodies confined within 120 AU and blown 



out by radiation pressure (Lecavelier des Etangs et al 



1996) 



- '/ ' 


/3 = 0.450 : 

(3 P , = 0.405 : 

(3 pr = 0.360 " 

/3 pr = 0.270 ■ 


: \ if/' i ~~ — — - 

<m i i . 






1 The Maxwell-Garnett effective medium theory is used to 
compute the complex index of refraction of an aggregate. 



100 200 300 400 

Distonce from the star [AU] 

Fig. 2. Surface density of the PB disk (bold-dashed line) 
that gives a good fit to the observed surface brightness 
distribution in scattered light. The resulting dust surface 
density integrated over the grain size (or /? pr ) distribu- 
tion is displayed with a bold-solid line. Both profiles are 
normalized to 100 at their maximum. The relative grain 
surface densities (thin lines) for five values of /3 pr ranging 
between 0.135 and 0.45 are also superimposed for compar- 



Applying the observational constraints described 
above, we tried several PB body surface density distri- 
butions. We were able to find a solution (Figure ||) that 
gives a fit to the scattered light images (Figure ||) . 

In this solution, the distribution of PB is confined in- 
side ~150 AU. Due to the radiation pressure acting on the 
grains, the outer part of the disk is filled with small par- 
ticles. This can be seen in Figure ||, where we have plot- 
ted the resulting surface density distributions for grains 
of different sizes. The final surface density of the grains 
averaged over the grain size distribution peaks at about 
110-120 AU from the star, close to the assumed outer edge 
of the parent disk (Figure ||) . 

At distances greater than about 30 AU, the main shape 
of the mid-plane surface brightness distribution in scat- 
tered light from Heap et al. (200C ) is well matched de- 
pending somewhat on the precise anisotropic scattering 
properties (Figure ||) . 

A radial power law index between —5 and —5.5 is 
measured at distances larger than 120 AU whe n the dust 
disk is seen almost edge-on as predicted by Lecavelier 
|des Etangs et al. (1996| ) and indeed observed with the 
HST/STIS instrument. 

Importantly, due to the effect of radiation pressure, the 
distribution of the grain sizes is a function of the distance 
from the star as shown in Figures | and |. 
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Fig. 5. The warped disk n ormalized for each vertical cut to its maximum surface brightness. To compare with Figures 
10 from Heap et al. (2000|) and to emphasize the vertical asymmetry, contours at 99%, 50% and 10% of the maximum 
surface brightness have been superimposed. 




40 60 80 100 200 

Distance from the star [AU] 

Fig. 3. Surface brightness distributions in mag.arcsec -2 
along the mid-plane of the disk in scattered light (R band) 
assuming the PB surface density shown in Figure ^. The 
solid line assumes isotropic scatte ring properties whereas 
the dash-dotted line assumes a Henyey Sz Greenstein 
(194l[j]phasc function with an asymmetry factor \g\ of 0.5 
The two dashed lines represent radial power laws with in- 
dices -1.1 and -5.5 as measured on the HST/STIS images 



3.3. Small scale vertical asymmetry (the warp) 

Due to the effect of the assumed perturbing planet, the 
precession of the orbital planes of PB within 70-80 AU 
results in a distribution with average inclination coinci- 



dent w lith the planetary orbit. Beyond 70-80 AU and up 
to the assumed outer edge of the PB disk (close to 140- 
150 AU), the PB orbits have undergone only partial pre- 
cession. Therefore, the mean inclination of their orbital 
planes to that of the initial disk ranges between 3° and 

about 0.5°. 

Mouillet et al. (1997b) showed that the warp, or the 
short scale asymmetry of the vertical position of the max- 
imum surface brightness, can be reproduced assuming all 
particles followed the distribution of the parent body disk. 
Adding the effects of radiation pressure actin g on grains 
with a giv en size distribution to the model of Mouillet et 
al. (1997b ) causes a modest shift of the warp position, and 
so does not affect the conclusions on the propagation and 
observability of the warp as seen in Figure ^|, which can 



30 AU 

1 30 AU 

1 50 AU 

180 AU 

210 AU 

250 AU 

300 AU 




150 200 250 

Distance from the star [AU] 

Fig. 4. Upper panel: Grain size distribution at different 
distances from the star. In each case the distribution is 
normalized such that the integral over a has been nor- 
malized to 1. Lower panel: N index assuming that the 
grain size distribution \dn(a)/da\ is proportional to: (1 — 
x ) x x~ N with x = a/a m [ n and a m i n ~ 0^=0.5 = 2K 
(see Table 0). In this example, M is between 1.4 and 2 
inside ~ 150 AU and up to 4 at larger distances. 



be directly compared to the HST/STIS images (Heap et 



al. 2000) 



3.4. The large scale vertical ( "butterfly" ) asymmetry 

We now investigate whether our model can account for 
the "butterfly" asymmetry. Basically, radiation pressure 
greatly extends the surface density distribution of the 
smallest particles, i.e. those which have the larger /3 pr val- 
ues, beyond that of the PB disk (figure [2]). Accordingly, the 
asymmetry in the vertical location of the surface bright- 
ness maximum originating in the PB disk can be trans- 
mitted outwards as illustrated in Figure ^. The result is a 
flared-like disk, with the so-called butterfly asymmetry. 

For illustrative purpose, we focus on the particular case 
for which the line of nodes of the planetary orbit and the 
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Fig. 7. Modeled disk at large distances from the star. A large scale asymmetry with respect to the mid-plane is found, 
similar to the observed "butterfly" asymmetry. 
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-600 -400 -200 200 400 600 

Distance from the star [AU] 

Fig. 6. The location of the maximum surface brightness 
for a vertical cut as a function of distance, for disks com- 
posed of grains with a single size as indicated, seen in scat- 
tered light. The line for /? pr = 0.001 traces the PB disk. 
The line for /3 pr = 0.45 indicates the way the small scale 
asymmetry (warp) can be extended outwards as a result of 
radiation pressure. The position of the observed mid-plane 
to the original mid-plane of the PB disk is superimposed 
in this particular case (line of nodes of the planet aligned 
with the line of sight and disk seen edge-on) . 



on the inclination of the planet to the PB disk mid-plane 
(3° in this model), and on the line of nodes of the planet 
to the line of sight. 



200 300 400 

Distance from the star [AU] 



Fig. 8. Ratio of half widths at half maximum above and 
below the observed mid-plane for a vertical cut as a func- 
tion of distance and for the extension of the disk seen in 
Figure 0. 



line of sight are aligned (the observability of the warp as 
a function of the location of the line of nodes is discussed 
in iMouillct ct al. 1997b| ). The simulated disk in scattered 



light then convincingly exhibits a large scale asymmetry 
as observed in STIS data up to 200 AU (Figure |). 

It is worth noticing that this asymmetry is present 
much further out than 200 A U (Figure |7|), in agreem ent 
with the large scale images of [Kalas fc Jewitt (1995] ) for 
instance. Figure || quantifies the radial dependence of the 
asymmetry by the ratio of half widths at half maximum 
above and below the observed mid-plane. The ratio evi- 
dences the position of the warp at 70-80 AU and shows 
that the disk remains asymmetric up to at least 500 AU 
as observed. In other words, the asymmetric distribution 
of PB responsible for the observed warp at 70-80 AU, and 
assumed to result from the perturbation of an inclined 
planet, leads naturally to a large scale asymmetry very 
similar to the "butterfly" asymmetry. Therefore, the ob- 
servability of the "butterfly" asymmetry depends as well 



3.5. Summary 

The dynamical approach we have followed allows us to 
reproduce in a consistent way most observations in scat- 
tered light, namely: the surface brightness distribution, 
the small scale (warp) and large scale ( "butterfly" ) asym- 
metries, assuming a planetary perturber and taking into 
account the effects of radiation pressure acting on the 
grains. But as expected, the model does not explain the 
NE-SW asymmetries. 

4. Thermal emission data 

4.1. IRAS fluxes 

A fully consistent disk model should reproduce both the 
scattered light images and thermal data. Our aim is to 
estimate the thermal IR (12, 25, 60 and 100 /jm) fluxes 
produced by the distribution found in Section ||. Table || 
summarizes the computed fluxes assuming the dust model 
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Table 1. Thermal emission in IRAS bands for different dust models. According to Li & Greenberg (1998), dust models 
with P — 0.98 and P — 0.95 assuming icy grains give a satisfactory fit to the overall SED. The percentages of ice have 
been adjusted to respect the ice to silicate core-organic refractory mantle mass ratio of 1 from Li & Greenberg (1998). 
Dust grains with P = 0.5 have a mass ratio of 0.25. Dust models with non-icy grains are also given for comparison. In 
each case, the dust surface density has been normalized in order to reproduce the scattered light profile and then the 
integrated thermal emission computed. With this procedure, the computed IR fluxes then also depend indirectly on the 
grain anisotropic scattering properties. There, we assume isotropic scattering properties, but the fluxes are increased 
by a factor of about 1.4 if the asymmetry factor \g\ increases to 0.5 and the ratio of the observed to computed flux is 
accordingly decreased by the same amount. 





IRAS 


P = 0.98 
0% ice 


P = 0.98 

4% ice 


P = 0.95 

0% ice 


P = 0.95 
10% ice 


P = 0.5 

0% ice 


P = 0.5 
50% ice 


Flux at 12 /im [Jy] 
Flux at 25 /im [Jy] 
Flux at 60 /im [Jy] 
Flux at 100 ^m [Jy] 


1.64±0.1 
10.1±0.5 
18.8±1 
11.2±1 


0.040 
3.5 
12.5 
10.4 


0.056 
2.1 
16.3 
10.2 


0.037 
3.4 
12.4 
10.3 


0.047 
1.9 
16.3 
10.2 


0.017 
2.4 
11.9 
9.8 


0.017 
1.6 
13.6 
9.7 


0,^=0.5 = 2K H | ™ 1Ce 




92.7 


92.7 
46.5 


37.1 


37.1 
18.9 


3.7 


3.7 
3.0 


Ratio of observed to 
computed flux at 12 /im 
further than ~30 AU (NE side) 
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proposed by Li & Greenberg (1998) for different porosities 



and amounts of ice. 

The 100 /im and 60 /im emissions (the latter slightly 
depending on the amount of ice) are correctly matched by 
the dust ring made of grains at the blow-out size limit that 
reproduces the scattered light profile. Such a consistency 
strongly supports our dynamical approach. 

On the other hand, the simulated 12 /im and 25 /im 
fluxes are smaller than observed. This shows that an ad- 
ditional population of hotter grains is required, in total 
agreement with previous works showing that most of the 
12 /im emission comes from a region close to the star (see 
for example Knacke et al. 1993), which is also con firmed 
by the resolved images from Pantin et al. (1997 ). This 



close population also certainly contributes to the inte- 
grated 25 /im flux. 



4.2. The 12 [im resolved images 



Most of the 12 /im flux originates within ~20 AU from the 



star as seen before. However, as shown by Pantin et al 



(1997 ), a small part of the 12 /im flux comes from regions 
up to ^100 AU. Our aim is here to see whether our model 
accounts for this "outer" 12 /im emission. 

Simulated and observed 12 /im images compare well 
in terms of radial shape down to about 20-30 AU (and 
especially when the most porous grains are not too icy, 
Figure ^|) but not in terms of flux. The ratio of observed 
to computed flux is given in the last line of the Table [l] for 
different dust properties. This shows that an additional 
population of efficient 12 /im emitters, and at the same 
time bad scatterers at near IR wavelengths, is needed be- 
tween 20-30 AU and 100 AU. Very small, typically submi- 
cronic, grains are a priori good candidates. For example, 




40 60 80 100 

Distance from the star [AU] 

Fig. 9. Linear brightness distribution in thermal emission 
at A = 12 /im assuming the PB surface density distribu- 
tion shown in Figure || and non-icy grains larger than 2K 
with P = 0.95. At each distance from the star, this linear 
brightness has been computed by multiplying the maxi- 
mum of flux along the vertical axis by the vertical thick- 
ness (in arcsec) of the disk at 12 /im given by the Figure 
4 from Pantin et al. (1997 ) and available for distances 



smaller than 70 AU. To compare with the Pantin et al 



(1997) resolved emission, the predicted fluxes have also 



been multiplied by a factor of 7, indicating that most of 
the grains that contribute to the 12 /im emission at dis- 
tances larger than 20-30 AU are not taken into account. 



the emission/absorption at 12 /im to scattering at 1 /im 
efficiency ratio is a few hundred for 0.1 /im grains with 
porosity P = 0.95 and 3 orders of magnitude larger if the 
: /im in size. 



grains are only 10 2 



4.3. Possible scenarios to account for the 12fim 
missing flux at large distances 

Core-mantle grains smaller than a/3 pr= o.5 = 2K have /3 pr 
ratios larger than 0.5. A close inspection of Figure || indi- 
cates that particles with f] pl ratios larger than typically 0.4 
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have quite similar surface densities and tend to mimic that 
of the /3 pr = 0.45 grains which are almost just unbound. 
In the following, we will use this property to simulate a 
dust disk with very small particles. 

We propose hereafter possible scenarios a priori able 
to explain the presence of very small grains at large dis- 
tances, and we check whether or not they may quantita- 
tively reproduce the 12 /im flux at large distances under 
reasonable assumptions. 

4.3.1. First scenario: Small grains produced by 
cometary evaporation 

A first possibility is that the submicronic grains produced 
close to the star, for instance by comet evaporation, re- 
sponsible for most of the integrated 12 /im flux and re- 
quired to explain the 10 /*m silicate spectroscopic features, 
also supply the disk beyond — 30AU with enough small 
dust particles through radiation pressure. This would give 
an attractive solution to the lack of 12 /im emission at 
distances larger than —30 AU in our model. 

We therefore assumed a population of very small grains 
close to the star and let them evolve under radiation pres- 
sure. It appears that, whatever the assumptions on the 
surface density distribution of the inner population, the 
radial shape of the 12 /im flux distribution does not match 
the one found by jPantinct al. (1997| ) beyond -30 AU. This 
is mainly due to the radial decrease of the grain temper- 
ature with distance from the star which tends to result 
in most of the emission at 12 /im coming from close to 
the star leading cither to very to small contributions fur- 
ther than —30 AU or to radial profiles for the emission not 
consistent with the resolved observations. 

A population of submicronic grains certainly exists 
close to the star but our modeling indicates that additional 
grains produced further than —30 AU should be present to 
account for the radial shape of the 12 /im emission. 

4.3.2. Second scenario: Small unbound grains below 
the blow-out limit 

We will now consider the dust disk that reproduces the 
scattered light data and extrapolate the initial a~ 35 grain 
size distribution to grains with f3 pT ratios larger than 0.5. 
We consider non-icy grains for simplicity and P = 0.95 
which gives good fits to the large wavelength fluxes. 

In other terms, we simply extend the grain size distri- 
bution that successfully fitted the scattered light images in 
Section^ to a minimum grain size of about 0.1 /im instead 
of a few tens of micrometers. This actually allows us to 
fit the 12 /im emission further than — 20 AU. The 60 /im 
and 100 /im integrated fluxes are unchanged as well as the 
scattered light profile inside 120 AU because the amount of 
large grains is not affected by the reduction in minimum 
grain size. However, beyond 120 AU, the scattered light 
surface brightness profile changes: it follows now a r~ 35 
radial power law (instead of r~ 5 ). r~ 3 - 5 is not consistent 



with STIS data, but is in better agreement with previous 



data (e.g. 


Golimowski ct al. 1993; 


Kalas & Jcwitt 1995 


Mouillct ct al. 1997b 


)■ 



It is therefore premature to conclude whether or not 
this extrapolation is correct on the sole basis of the sur- 
face brightness distribution at large distances. Moreover, 
the rate of production of such small grains must be con- 
tinuous and more efficient than indicated by a simple ex- 
trapolation from the larger grains, since being unbound, 
their life-time in the disk is very short, in fact 1 to 10 3 
times shorter than grains with sizes just above the blow 
out limit. The mass loss rate from the system due to the 
blow out of these grains is also correspondingly larger. 

4.3.3. Third scenario: Small compact bound grains 
produced by collisions 

The physics of aggregates, present in circumstellar 
environments, is complex and not well understood. 
Catastrophic disruption of an aggregate can occur if the 
impact velocity is high enough (ty pically a few tenths of 
km.s -1 , Dominik and Ticlens 1997, for instance), produc- 
ing nearly compact submicronic grains. In the vicinity of 
(3 Pictoris, such grains smaller than — 10 -2 /im may have 
/3 pr ratios smaller than 0.5 depending on their precise 
chemical composition, i.e. if they are not too refractory 
(e.g. Artymowicz 1988| ). This leads ultimately to a discon- 
tinuity in the grain size distribution which could appear 
depleted in grains of sizes ranging between — 10~ 2 /im and 
2K. 

We tested this scenario by assuming an additional pop- 
ulation of 10~ 2 /im compact silicate grains with a sur- 
face density distribution that corresponds to that of the 
/3 pr = 0.45 particles. The 12 /im resolved emission is very 
well matched in terms of shape and flux —20-30 AU with a 
mass of 10~ 2 /im grains of the same order of magnitude as 
the mass contained in the particles of size 2K. The far-IR 
emission and the scattered light images, in particular the 
r~ 5 radial shape of the surface brightness distribution fur- 
ther than 120 AU, are not affected by the presence of such 
very small grains since they are very poor scatterers. The 
presence of such very small grains, and if so, their chem- 
ical composition as well as their size distribution, remain 
open issues. 



5. Summary and conclusions 



As in Mouillet et al. (1997b), we assumed a disk of plan- 
etesimals orbiting around (3 Pictoris and subject to the 
gravitational perturbation of a planet on an inclined or- 
bit. 

Assuming that collisions among these planetesimals 
produce grains with a size distribution proportional to 
a -3.5 between the blow out limit and millimeter sizes, and 
assuming that these grains are subject to radiation pres- 
sure in addition to gravitational forces, and using plausible 
chemical composition and optical properties for the grains, 
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we are able to reproduce the scattered light images (inten- 
sity, radial distribution, warp, "butterfly" asymmetry) as 
well as the far IR integrated fluxes. The PB are located be- 
tween 20 and 150 AU with a peak at -100-120 AU. Due to 
radiation pressure, the grain radial distribution strongly 
depends on their size and the resulting dust disk extends 
very far from the star. The "butterfly" asymmetry is noth- 
ing more than the mirror of the warp of PB disk at larger 
distances due to radiation pressure. 

The 12 /im images require an additional population 
of small grains close to the star as was evident several 
years ago and attributed to comet evaporation. We have 
shown that radiation pressure acting on these grains 
does not produce a distribution that can account for the 
observed radial sha pe of the 12 /im emission as found by 
Pantin et al. (1997 ) at distances larger than —20-30 AU. 
A population of grains smaller than (or comparable 
to) the blow-out limit in addition to the population 
necessary to reproduce the scattered light images may 
account for this > 30 AU 12 /im emission. It is found 
that the extrapolation of the a -3 5 distribution down to 
such small sizes is somewhat problematic. A tentative 
alternative could be that there is a discrete population 
of grains with size less than 10~ 2 /im in addition to 
the main population of grains. This is plausible as we 
assume that the main population of grains is made 
of porous aggregates ultimately breaking down and 
producing < 10 -2 /im grains that could be responsible 
for the emission further than 30 AU at 12 /zm, without 
affecting the scattered light images and the 60-100 /im 
integrated fluxes. As another alternative, it is also 
possible to argue that the a~ 3 5 law is unlikely close to 
the blow-out limit: grains close to the blow-out size may 
be over-abundant simply because of the lack of smaller 
grains able to efficiently destroy them. In this latter 
scenario, a larger number of small grains is explained as a 
result of collisions, without involving additional processes. 

Finally, our consistent modeling of the (3 Pictoris disk, 
combining the various available observations has the fol- 
lowing limitations: 1/ the solution found is probably not 
unique as various parameters (optical properties, size dis- 
tribution, spatial distribution, grain production mech- 
anisms) are intricately involved; 2/ the application of 



the 



-3.5 



size distribution all the way between blow-out 



size and parent bodies is questionable; 3/ the side- to- 
side asymmetry in the 12 /im and 20 /im images has not 
yet been explained. Obviously, this asymmetry originates 
close to the star and might be due to the evaporating 



comets as was proposed by Lecavelier des Etangs et al 
(19960 

Further investigations will be able to benefit from new 
direct constraints on the spatial distribution of grains, as 
observed at various wavelengths (scattered light and ther- 
mal emission high angular resolution images), for the inner 
disk (< 50 AU), and optical properties of the grains. 
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